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Abstract:  The  electronic  band  structures  and  optical  properties  of  type-II 
superlattice  (T2SL)  photodetectors  in  the  mid-infrared  (IR)  range  are  in¬ 
vestigated.  We  formulate  a  rigorous  band  structure  model  using  the  8-band 
k  ■  p  method  to  include  the  conduction  and  valence  band  mixing.  After 
solving  the  8  x  8  Hamiltonian  and  deriving  explicitly  the  new  momentum 
matrix  elements  in  terms  of  envelope  functions,  optical  transition  rates  are 
obtained  through  the  Fermi’s  golden  rule  under  various  doping  and  injection 
conditions.  Optical  measurements  on  T2SL  photodetectors  are  compared 
with  our  model  and  show  good  agreement.  Our  modeling  results  of  quantum 
structures  connect  directly  to  the  device-level  design  and  simulation.  The 
predicted  doping  effect  is  readily  applicable  to  the  optimization  of  photode¬ 
tectors.  We  further  include  interfacial  (IF)  layers  to  study  the  significance  of 
their  effect.  Optical  properties  of  T2SLs  are  expected  to  have  a  large  tunable 
range  by  controlling  the  thickness  and  material  composition  of  the  IF  layers. 
Our  model  provides  an  efficient  tool  for  the  designs  of  novel  photodetectors. 
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1.  Introduction 

Antimony-based  type-II  superlattices  (T2SLs)  have  attracted  great  interest  from  researchers 
due  to  their  wide  range  of  applications  in  optoelectronic  devices  for  the  past  decades.  Such 
periodic  quantum  structures  were  first  proposed  by  L.  Esaki  [1]  in  1977  and  applied  to  in¬ 
frared  (IR)  applications  by  Smith  and  Mailhiot  [2]  in  1987.  The  broken-gap  band  alignment  in 
type-II  heterojunctions  allows  a  narrow  effective  band  gap  between  the  quantized  conduction 
band  (CB)  and  valence  band  (VB)  states.  Long  optical  transition  wavelength  up  to  mid-IR  (3 
jlm  to  30  flm )  becomes  achievable  for  interband  transitions  without  resorting  to  intersubband 
transitions  [3].  The  type-II  material  system,  a  promising  candidate  for  active  layer  designs  in 
mid-IR  optoelectronic  devices,  has  recently  demonstrated  room  temperature  operation  in  both 
interband  cascade  lasers  [4]  and  photodetectors  [5].  Antimony -based  III-V  compound  SLs  have 
higher  growth  quality  and  uniformity  than  II-VI  compounds  such  as  HgCdTe  alloys,  and  a  large 
reduction  of  defects  is  expected.  A  large  tunable  range  of  operation  wavelength  can  be  obtained 
by  varying  layer  thicknesses  rather  than  precise  control  of  mole  compositions.  Suppression  of 
Auger  recombination  is  also  predicted  in  antimony -based  T2SLs  [6]. 
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The  k  •  p  method  [7-10]  serves  as  an  efficient  model,  which  covers  from  material  level, 
quantum  structural  level,  up  to  the  device  level  modeling.  Despite  the  loss  of  atomistic  de¬ 
tails  [11],  this  method  is  desirable  for  large-scale  problems,  long  wavelength  applications,  and 
device  modeling  [9, 12, 13].  The  band  structure  near  high-symmetry  extremum  points  (e.g.  T,  X, 
L)  is  solved  through  the  perturbation  theory,  and  the  output  information  is  accurate  enough  for 
simulating  optoelectronic  processes  near  the  semiconductor  band  gap.  The  results  from  the  k  •  p 
method  can  be  readily  used  for  device-level  analysis  and  design.  Our  previous  work  has  applied 
this  method  for  the  gain  optimization  in  type-II  quantum  cascade  lasers  [9].  The  calculated  ma¬ 
terial  gain  spectra,  which  enter  into  the  rate  equations  for  light-emitting  devices,  allow  for  more 
accurate  modeling  of  the  light  output  vs.  injection  current  (L-I)  curves.  We  have  also  analyzed 
the  T2SL  photodetector  performance  [13]  including  the  quantum  efficiency  and  responsivity, 
based  on  an  analytical  drift-diffusion  model  with  the  input  of  the  calculated  absorption  spectra. 

In  this  paper,  we  formulate  a  rigorous  model  for  band  structures  and  optical  properties  using 
the  8-band  k  •  p  method.  We  derive  the  new  8-band  momentum  matrix  elements  for  TE  and  TM 
polarizations.  Our  expressions,  explicitly  derived  in  terms  of  the  envelope  functions,  are  more 
straightforward  than  those  by  Szmulowicz  [14]  (which  were  expressed  in  terms  of  matrices  of 
operators),  and  more  general  than  those  by  Chang  and  James  [15]  (which  were  derived  only  for 
intersubband  transitions  in  the  valence  band,  ignoring  the  CB-VB  mixing).  Our  results  under 
the  axial  approximation  [9]  significantly  simplify  the  computation  and  are  readily  applicable 
for  modeling  optical  transitions  with  CB  and  VB  mixing,  while  the  ab-initio  band  structure 
models  are  not  yet  applied  for  modeling  T2SL  optical  properties  and  device  performance. 

One  challenge  in  the  k  ■  p  modeling  is  the  determination  of  the  input  bulk  band  edge  pa¬ 
rameters.  The  electron  effective  mass  and  the  Luttinger  parameters,  which  account  for  the  hole 
effective  masses,  become  sensitive  parameters  when  modeling  narrow  gap  materials,  such  as 
InAs  and  InSb  [16].  Band  edge  parameters,  which  enter  into  the  Hamiltonian  formulation,  need 
to  be  corrected  from  experimental  values.  Such  correction  is  only  necessary  for  the  8-band  k  •  p 
method  or  beyond,  since  the  CB-VB  coupling  is  already  included  in  the  off-diagonal  terms  in 
the  8x8  Hamiltonian.  The  Luttinger  parameters  are  also  corrected  due  to  the  change  of  the 
definition  for  class-B  states  (remote  bands)  in  the  Lowdin’s  perturbation  method  [17].  In  the  6- 
band  case  [18],  s-like  CB  states  are  considered  as  “remote  bands”,  while  in  the  8-band  case  [7], 
s-like  CB  states  belong  to  the  8-fold  bases  in  class-A  states.  We  present  the  correction  formulas 
and  summarize  the  corrected  parameters  for  commonly  used  type-II  materials. 

We  verify  our  model  using  two  different  boundary  conditions:  Dirichlet  and  Periodic  Bound¬ 
ary  Conditions  (DBC  and  PBC).  Important  optical  properties  such  as  absorption  coefficient  and 
spontaneous  emission  rate  are  calculated  and  compared  with  the  experimental  data.  Our  theoret¬ 
ical  results  on  absorption  and  photoluminescence  (PL)  spectra  agree  very  well  with  the  optical 
measurements.  The  measured  quantum  efficiency  (QE)  spectra  of  a  mid-wave  and  a  long-wave 
infrared  photodetectors  [19,20]  serve  as  verifications  of  our  model.  We  further  demonstrate 
the  doping-dependence  of  the  T2SL  optical  properties.  The  light  p-doping  does  not  change  the 
cutoff  wavelength  and  slightly  affects  the  quantum  efficiency.  This  allows  us  the  freedom  to 
optimize  T2SL  detectors  using  different  levels  of  p-doping. 

Another  challenge  of  modeling  heterostructures  is  the  determination  of  the  band  alignment, 
or  the  commonly  used  parameter:  valence  band  offset  (VBO).  Variation  in  the  measured  VBO 
causes  a  large  fluctuation  in  the  modeled  results  for  narrow-gap  materials.  The  interfacial  (IP) 
layers  formed  at  the  heterojunctions  play  an  important  role  in  the  uncertainty  of  band  align¬ 
ments  besides  experimental  errors.  Previous  modeling  work  mostly  included  the  IP  effect  by 
adjusting  the  VBO  values  or  using  a  graded  potential  profile.  However,  these  approaches  failed 
to  account  for  the  large  IP  lattice  mismatch  and  the  strain  effect  on  the  band  structure  and  the 
optical  properties.  In  this  paper,  we  include  actual  IF  layers  in  our  model  and  demonstrate  the 
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significance  of  the  IF  effect.  IF  layers  are  often  used  to  study  lattice  mismatch  and  achieve 
strain  balance  [21-24].  Beyond  previous  studies,  we  predict  with  our  model  that  one  can  obtain 
a  large  tunable  range  of  optical  properties  by  controlling  the  IF  layer  type  and  thickness. 


2.  Theoretical  principles 

2.1.  The  8-band  k  p  method  with  the  effective  mass  correction 

In  this  paper,  the  band  structure  and  optical  properties  of  T2SLs  are  modeled  with  the  8-band 
k  ■  p  method,  which  is  based  on  the  Kane’s  formulation  [7, 10],  the  Luttinger-Kohn  (LK)  model 
[8,  10]  for  degenerate  bands,  and  the  Pikus-Bir  theory  [10,25]  for  strain  effects.  Due  to  the 
narrow  band  gap  of  type-II  materials,  it  is  necessary  to  use  an  8-band  coupled  Hamiltonian  to 
include  the  strong  interaction  among  the  8  bands  nearest  to  the  band  gap,  namely  the  conduction 
band  (CB),  the  heavy-hole  (HH),  the  light-hole  (LH),  and  the  spin-orbit  split-off  (SO)  bands, 
each  with  a  double  spin  degeneracy.  The  8x8  LK  Hamiltonian  built  under  these  8  bases  are 
first  block-diagonalized  into  two  decoupled  4x4  upper  and  lower  Hamiltonian  matrices  as  [9] 


ijLK  _ 

^8x8  — 


^4x4  (^f) 


H^t) 


(1) 


using  the  basis  transformation  from  Eq.  (A-3)  to  Eq.  (A-4).  The  band  structure  is  generally 
anisotropic  in  the  transverse  plane  and  depends  on  the  direction  of  the  transverse  wave  number 
kf.  After  the  axial  approximation,  the  Hamiltonian  matrix  elements  and  the  envelope  functions 
depend  only  on  the  magnitude  of  kf .  The  upper  and  lower  4x4  Hamiltonians  are  written  explic¬ 
itly  in  Eq.  (A-l).  The  upper  and  lower  eigen-states  of  the  m-  and  n-th  subbands  can  be  written 
as  a  linear  conbination  of  the  new  basis  functions  modulated  by  the  envelope  functions  as, 


^(k„r) 

^(k„r) 


t-Tt  4 

—TffihSm \kt,z))\ui), 

v  A  i=i 
pte-t-r,  8 

-JffXSn  {kt,z))\ui). 

vA  ,=5 


(2) 


The  eigen-energies  and  eigen-states  at  a  given  kt  can  then  be  solved  from  the  matrix  equation. 


’  g{1Hkt,z)  ' 

'  gW(kt,z)  ' 

g(2\k„z) 

=  Eu  (kt) 

g{2\kt,z) 

g(3)(£r,z) 

g(3>(k,,z) 

.  g[A\ktlz)  _ 

.  g(4)(h,z )  . 

for  the  upper  Hamiltonian.  Similar  equation  holds  for  the  lower  Hamiltonian  H^y4(kt).  At  the 
zone  center  ( kt  =  0),  there  is  usually  a  dominant  component  among  the  envelope  functions 
which  the  band  can  be  named  after.  For  wide  band  gap  materials,  the  CBs  are  weakly  coupled 
to  the  valence  bands  (VBs)  and  usually  solved  separately  from  the  6x6  Hamiltonian  with 
VB  mixing  [18],  and  the  weak  CB-VB  coupling  is  included  in  the  effecitve  mass.  In  the  8x8 
Hamiltonian  for  narrow  band  gap  materials,  the  electron  effective  mass  in  the  diagonal  matrix 
elements  need  to  be  corrected  since  the  off-diagonal  elements  already  include  the  coupling.  The 


#156797  -  $15.00  USD  Received  19  Oct  2011;  revised  29  Nov  2011;  accepted  19  Dec  2011;  published  18  Jan  2012 
(C)  2012  OSA  16  January  2012  /  Vol.  20,  No.  2  /  OPTICS  EXPRESS  2322 


corrected  effective  mass,  denoted  by  m'c,  enters  into  the  eight  diagonal  matrix  elements  as, 


(Mi,5|^8x8lMi,5)  =EC+Ae  +  f  +A1  J  (kf  +k\)  =  Ec  +Ae  +  ^y(k2  +^z)i 

Tr 

(m2,6|^8xsI«2,6)  =EV-Pe-Qe~  —  [(yi+y2)^  +(7i  -2y2)*f] , 

h 2 

(M3, 7 |^8x8 1^3,7)  =EV-Pe+Qe-^-  [(71-72)^?  +(7l+272)^?]  > 

h2 

{u4$\H%x8\uw)  =EV-Pe-  A-  —71  (k2  +  k2) . 


(4) 


where  Ec  and  Ev  are  the  unstrained  bulk  CB  and  VB  edges,  and  A  is  the  spin-orbit  split-off 
energy.  Ae,  PE,  and  Q£  are  the  strain  contributions,  which  are  expressed  explicitly  in  Eq.  (A-2). 

The  corrected  effective  mass  can  be  written  in  terms  of  Kane’s  paramemter  A',  which  in¬ 
cludes  the  momentum  matrix  elements  between  the  s-like  CBs  and  remote  bands  with  T^  sym¬ 
metry  type  [7] . 


A'  = 


ft2  y  \(S\p.x\uj)\2 
mlj  Ec-Ej 


(5) 


The  original  effective  mass  mc,  which  can  be  obtained  from  experiments,  describes  the  curva¬ 
ture  of  the  eigen-energy  dispersion.  The  expression  [26, 27]  for  mc  can  be  derived  by  solving 
energy  dispersion  E(kt)  near  the  zone  center  with  expansion  up  to  0(k2), 


mo  =  2mp  ,  Eg  +  2A/3 

mc  h2  PEs(Eg  +  A)' 


(6) 


Therefore  the  correction  of  electron  effective  mass  from  the  experimental  values  is  written  as, 


mo  mo  Eg+2A/3 

m!c  mc  p  Eg  ( Eg  +  A) 


The  Luttinger  parameters  obtained  from  experiments  correspond  to  the  Dresselhaus  [28]  defi¬ 
nition,  where  the  class-A  states  for  Lowdin’s  perturbation  method  [17]  use  the  6-fold  VB  bases. 
However,  for  the  8-band  LK  Hamiltonian,  the  class-A  states  should  be  defined  using  the  8-fold 
bases  [7,26,27]  including  the  two  s-like  CBs.  Therefore  the  extra  coupling  term  between  the 
s-like  CB  states  and  the  p-like  VB  states  needs  to  be  subtracted  from  the  original  Luttinger 
parameters  [7,26-29], 


The  original  and  the  corrected  material  parameters  are  shown  in  Table  1.  Values  are  updated 
from  [9]  to  include  the  corrections  and  the  original  parameters  are  within  error  margins  in  [29]. 
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Table  1 .  Input  material  parameters  for  the  8-band  k  •  p  method 
Parameter  Symbol  (unit)  GaSb  InAs  InSb  AlSb 


Lattice  constant  (A) 

a 

6.0854 

6.0522 

6.4717 

6.1297 

Deformation  potentials  (eV) 

ac 

-7.5 

-5.08 

-6.94 

-4.5 

av 

0.8 

1 

0.36 

1.4 

b 

-2 

-1.8 

-2 

-1.35 

Elastic  stiffness  constant 
(10]1dyne/cm2) 

Cn 

8.842 

8.329 

6.847 

8.769 

Elastic  stiffness  constant 
(101  'dyne/cm2) 

Cn 

4.026 

4.526 

3.735 

4.341 

Band  gap  at  OK  (eV) 

Eg 

0.81 

0.42 

0.235 

2.386 

Band  gap  at  77K  (eV) 

Eg 

0.80 

0.407 

0.227 

2.374 

Electron  effective  mass 

mc/mo 

0.0412 

0.0224 

0.0135 

0.13 

Corrected  electron  effective  mass 

m'c/mo 

1.053 

2.51 

1.84 

2.716 

Luttinger  parameters 

rf 

11.84 

19.4 

32.4 

4.15 

4.25 

8.545 

13.3 

1.28 

73 

5.01 

9.17 

15.15 

1.75 

Corrected  Luttinger  parameters 

7\ 

4.86 

6.30 

16.78 

1.76 

72 

0.76 

1.99 

5.61 

0.085 

73 

1.52 

2.62 

7.46 

0.55 

Optical  matrix  parameter  (eV) 

Ep 

22.4 

21.5 

23.3 

18.7 

Spin-orbit  split-off 
energy  (eV) 

A 

0.81 

0.38 

0.81 

0.65 

Valence-band  offset 
energy  (w.r.t.  GaSb)  (eV) 

VBO 

0 

-0.56±0.01 

0.03 

-0.38 

2.2.  Dirichlet  boundary  condition  (DBC)  vs.  periodic  boundary  condition  (PBC) 

Let  the  thickness  of  one  superlattice  period  be  Lp  and  the  total  thickness  with  N  periods  be 
Lj  =  NLp.  With  DBC,  the  solution  region  is  Lj  and  the  wave  functions  at  both  ends  are  set  to 
zero,  or  equivalently  introducing  infinite  potential  barriers  outside  the  solution  region. 

g£\kt,-LT/2)=g(£{kt,LT/2)=0  (9) 

With  PBC,  the  solution  region  is  Lp  and  the  wave  function  values  at  both  ends  are  related  by 
Bloch’s  theorem.  An  extra  quantum  number  q  is  introduced  as  a  Bloch  wave  number  qz.  to 
account  for  the  superlattice  periodicity  in  z-direction.  The  PBC  is  formulated  as, 

f  gm  ( kt,z  +  Lp)=e‘qLpg<in\k,,z ) ,  ( l ^ 

\  gm(kt,z+LT)  =  g$(kt,z+NLp)  =  gm\kt,z). 

The  eigen-equation  is  obtained  as, 

eiNqLp  =  1,  and  q=—.  (11) 

H  NLp 
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Fig.  1.  Wavefunctions  at  k,  =  0  solved  using  the  Dirichlet  boundary  condition  (DBC)  and 
the  periodic  boundary  condition  (PBC).  Ig^1*]  of  the  lowest  Cl  band  (green)  and  |g^21|  of 
the  highest  HH1  band  (red)  are  plotted  with  Cl  and  HH1  band  edges  (blue  dashed).  The 
zero  energy  reference  is  chosen  as  the  unstrained  conduction  band  edge  of  InAs. 


Fig.  2.  Energy  dispersion  relation  to  transverse  wave  number.  Blue  dashed  lines  indicate 
the  band  edge  energies.  The  zero  energy  reference  is  chosen  as  the  unstrained  conduction 
band  edge  of  InAs. 


where  q  takes  N  (or  /V  +  1 )  values  in  the  first  Brillouin  zone  (—n/Lp,n/Lp), 


n 


0,  ±1,  ±2, ...,  ±—  N  is  even, 
0,  ±1,  ±2,...,  ±^~^~  N  is  odd. 


(12) 


The  solution  region  of  Lp  is  solved  for  N  (or  N  +1)  times  with  different  Bloch  phase  shift 
determined  by  q.  The  number  of  unknowns  for  PBC  is  N  (or  N  +  1)  times  fewer  than  that  of 
DBC  but  the  eigen-problem  is  solved  for  N  (or  N  +  1)  times.  Since  most  eigen-solvers  have  a 
complexity  greater  than  0{N2),  the  computational  cost  using  PBC  is  much  less  than  DBC  when 
N  becomes  large. 

After  choosing  the  boundary  condition,  the  coupled  eigen-equations  in  Eq.  (3)  are  solved 
numerically  using  the  finite-difference  method.  Using  DBC,  each  subband  (e.g.  Cl  or  HH1)  is 
splitted  into  N  coupled  states,  which  is  consistent  with  the  coupled-mode  theory.  The  band  edge 
energies  converge  as  N  increases.  Using  PBC,  the  band  edge  is  determined  by  the  q  —  0  states, 
where  wavefunctions  in  one  period  are  repeated  in  other  periods  without  any  phase  change. 

The  dominant  envelope  functions  (g' 1  for  CB,  g':2  for  HH)  are  plotted  in  Fig.  1.  The  zero 
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energy  reference  is  chosen  as  the  unstrained  conduction  band  edge  of  InAs.  Comparing  DBC 
and  PBC,  band  edge  energies  (blue  dashed)  agree  well.  The  thin  GaSb  layers  cause  strong  cou¬ 
pling  between  electrons  confined  in  the  InAs  layers.  For  DBC,  this  strong  coupling  brings  forth 
a  slow-varying  envelope  that  modulates  CB  wavefunctions,  which  is  mathematically  caused  by 
the  infinite  barriers  at  both  ends.  For  PBC,  the  strong  coupling  effect  is  reflected  in  the  large 
non-zero  mean  value  of  the  wavefunctions,  even  though  the  magnitude  of  the  wavefunctions 
repeats  for  all  periods.  Despite  the  difference  in  wavefunctions,  the  optical  transition  matri¬ 
ces  for  DBC  and  PBC  should  converge,  because  the  overlap  integration  can  be  approximately 
separated  into  a  product  of  two  integrals,  one  over  fast-varying  part  in  a  single  period  of  SL, 
and  the  other  over  the  slow-varying  part  across  the  whole  solution  region.  The  average  of  the 
slow-varying  envelope  is  equivalent  to  the  large  non-zero  mean  when  using  PBC. 

The  energy  dispersion  relations  E [kt)  are  shown  in  Fig.  2,  calculated  using  9  periods  of  SLs. 
The  assignment  of  the  HH  or  LH  band  names  is  according  to  the  dominant  envelope  function 
at  the  band  edge  (i.e.  A'  for  HH  and  g(  l)  for  LH).  The  effective  band  gap  and  the  band  edge 
energies  for  DBC  and  PBC  agree  well.  For  DBC,  each  subband  is  split  into  9  coupled  states, 
which  converge  at  higher  k, .  For  PBC,  each  subband  is  split  into  5  discrete  energy  levels  because 
of  the  double  degeneracy  among  the  9  different  q  values  (E+q(kt)  =  E~q(kt)). 


2.3.  Absorption  and  photoluminescence  spectra 

The  optical  absorption  coefficient  [10]  in  a  quantum  well  structure  can  be  derived  from  the 
Fermi’s  golden  rule  and  expressed  as. 


a{Tm)  =- 


ne~ 


y  y  _L  f2n  dA  [ 

ch,hLTJ o  2nJo 


nrce0m^(O 

x\fAm(kt)-f?l’n(kt)]L(kt,hco), 


’  ktdkt 
2k 


\(W^n\e-p\^m)\2 


where  the  Fermi  distribution  functions  are 


m) 


l+exp(£c’"(^ 


l+exp(£v’m(^)  Fv) 


(13) 


(14) 


for  electron  occupation  probability  in  the  «-th  conduction  subband  with  an  energy  Ecn(kt)  and 
the  ra-th  valence  subband  with  an  energy  Evm  (kt),  respectively.  Fc  is  the  quasi-Fermi  level  for 
the  conduction  band  electrons  and  Fv  is  the  quasi-Fermi  level  for  the  valence  band.  (1  —  f"‘ )  is 
the  occupation  probability  by  the  holes.  The  summation  over  G\  and  a?  includes  the  transitions 
among  eigen-states  of  both  upper  and  lower  4x4  Hamiltonian,  and  the  spin  degeneracy  is  also 
accounted  for. 

Evaluation  of  the  8x8  momentum  matrix  elements,  as  suggested  by  Chang  and  James  [15], 
and  Szmulowicz  [14],  is  different  from  that  of  the  6x6  case  [18],  where  the  mixing  among  VB 
states  is  included  but  not  the  mixing  between  the  CB  and  the  VB  states.  In  the  6x6  case,  all 
interband  transition  matrix  elements  are  proportional  to  (S\pz\Z)  since  the  CB  states  are  purely 
s-like  and  the  VB  states  are  mixtures  of  p-like  states.  In  the  8  x  8  case,  both  the  CB  and  the  VB 
states  include  the  8-band  mixing  and  interband  transition  is  allowed  from  p-like  to  p-like  states. 
Such  a  transition  (e.g.  (Z\pz\Z))  is  non-vanishing  when  Lowdin’s  perturbation  method  [17]  is 
applied  to  include  the  remote  bands  denoted  as  class-B.  Since  the  momentum  operator  does 
not  commute  with  the  Hamiltonian,  even  the  new  bases  {|  m,)}  cannot  simultaneous  block- 
diagonalize  the  momentum  matrix.  Such  a  dense  8x8  momentum  matrix  can  be  obtained  from 
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the  Hamiltonian  using  the  lsf-order  k  •  p  perturbation, 

P8x8  =  -^-Vk//8X8,  (15) 

since  «  Hq  +  (fi/mo)k  ■  p.  The  gradient  in  Eq.  (15)  is  expanded  by  the  2«oLorder  derivatives 
of  the  Hamiltonian.  Then  the  momentum  matrix  elements  can  be  written  as. 


(^‘•"(k<,z)|e-p|‘P?2’m(k f,z))=e-X  [ dz(g{j\kt,zj) 

i,j  J 


'"°V  r, 

LTVk//" 


gm\kt,z),  (16) 


which  is  consistent  with  [15]  and  [14].  We  further  derive  explicitly  the  8x8  momentum  matrix 
elements  for  TE  and  TM  polarizations  in  terms  of  the  envelope  functions  in  Appendix  B. 

The  direction  z  is  chosen  as  the  epitaxial  growth  direction  of  the  T2SLs.  Therefore  in  Eq. 
(16),  e_L z  and  e\\z  for  TE-  and  TM-polarized  optical  transitions,  respectively.  For  mid-IR  ab¬ 
sorption  spectra,  TE  polarization  is  of  main  interest  because  of  its  major  contribution  among 
top-incident  light,  especially  near  transition  band  edge,  where  normally  HH1-C1  is  dominating. 

The  lineshape  function  L{kuh(o)  uses  the  normalized  Gaussian  distribution  to  account  for 
a  finite  transition  linewidth  due  to  various  scattering  mechanisms,  with  y  being  the  linewidth 
parameter. 


L(ktfico) 


1 


(Ec,n(kt)  —EVtm(k,)  —Tl(0)2 
2  y2 


(17) 


Unlike  absorption  or  stimulated  emission  where  the  net  transition  rate  is  calculated,  spon¬ 
taneous  emission  considers  only  the  downward  transition  from  the  CBs  to  the  VBs.  It  only 
requires  the  initial  state  to  be  occupied  and  the  final  state  empty,  assured  by  the  term  fc(  1  — 

The  photoluminescence  (PL)  spectrum  is  proportional  to  the  spontaneous  emission  spectrum. 


nre2(0  V1  ^  [2K  d<t>  [°°  hdkt 

c37th£oml  0^“  LT  Jo  2k  Jo  2 n 

|  (V?  ’n\e-  p\^m)  1 2/eff|  ■"(*»)[!-  fva2’m(kt)]L(kt,  hco). 


(18) 


The  integrations  over  kt  in  Eq.  (13)  and  (18)  are  performed  numerically  by  truncating  at  a 
finite  range  (k,  =  0.2/A).  The  lineshape  function  approaches  zero  rapidly  when  k,  is  large, 
which  indicates  the  high-1:,  transitions  contribute  little  to  the  infrared  absorption  or  emission. 
For  absorption  in  an  intrinsic  T2SL,  the  quasi  Fermi-level  separation  is  zero  ( Fc  =  Fv)  and  the 
Fermi-level  lies  in  the  middle  of  effective  band  gap.  Therefore  the  Fermi  function  difference 
(fv  —  fc)  is  1 .  For  PL,  low-level  injection  creates  a  small  quasi  Fermi -level  separation  but  the  two 
quasi  Fermi-levels  are  both  deep  within  the  gap.  For  gain  spectrum,  the  population  inversion 
condition  requires  negative  Fermi  function  difference  (/,,  fc  <  0)  in  order  to  have  a  positive 
gain  value.  Therefore  the  injection  level  needs  to  be  high  enough  to  create  a  large  quasi  Fermi- 
level  separation  (Fc  —  Fv).  Gain  is  equal  to  negative  absorption. 


g(hco)  =  -a(hco). 


(19) 


For  both  PL  and  gain  modeling,  Fc  and  Fv  are  related  to  the  injected  carrier  density  by. 


n  = 


P  = 


(20) 
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To  solve  Fc  and  Fv  using  Eq.  (20),  we  obtain  carrier  densities  n  and  p  from  the  charge  neutrality 
condition. 


n+NA=p+N+,  (21) 

where  NA  is  the  ionized  acceptor  concentration,  Np  is  the  ionized  donor  concentration.  The 
carrier  densities  n  and  p  can  be  found  from  n  =  n o  +  8n,  p  =  po  +  Sp,  and  8n  =  8  p,  where 
no  and  po  are  the  electron  and  hole  concentrations,  respectively,  at  thermal  equilibrium.  8n  and 
8p  are  the  excess  electron  and  hole  concentrations  due  to  electrical  or  optical  injection.  For  PL 
or  low-level  injection  condition, 

„? 

no  =  —  <  8n  =  8p  <  po  =Na  -N%, 

Po 
n 2 

po  =  —  <C  8n  =  8p  <C  no  =  Np  —Na  , 
no 

no  =  PO  =  ni  <  Sn  =  Sp , 

For  absorption  spectra  at  zero  bias,  the  excess  carrier  injection  is  zero  ( 8n  =  Sp  =  0). 

3.  Simulation  results  and  discussion 

In  this  paper,  we  model  the  absorption  spectrum  and  the  PL  spectrum  for  an  intrinsic  T2SL  with 
an  InAs/GaSb  44  A/21  A  layer  structure.  Specific  interband  transitions  (e.g.  HH1-C1,  LH1-C1) 
are  calculated  together  with  the  total  transition  as  shown  in  the  inset  of  Fig.  3.  The  cutoff 
energies  for  different  types  of  interband  transitions  are  consistent  with  the  effective  band  gaps 
between  different  subbands  predicted  from  Fig.  2.  The  HH1-C1  transition  turns  on  at  around 
130  meV  in  Fig.  3.  The  HH1-C2  transition  appears  at  the  second  step  of  the  HHl-to-all  CBs 
absorption  curve  (black  curve  in  Fig.  3  inset),  which  agrees  with  the  effective  band  gap  between 
C2  and  HH1  subbands  (470  meV)  in  Fig.  2.  The  LH1-C1  starts  at  around  300  meV,  which  is 
the  effective  band  gap  between  Cl  and  LH1  subbands.  The  peak  of  LH1-C1  transition  appears 
at  500  meV  due  to  the  large  wavefunction  overlap  at  this  transition  energy  and  the  contribution 
from  higher  kt  states. 

The  external  quantum  efficiency  (QE)  of  an  InAs/GaSb  45  A/24  A  T2SL  photodetector  is 
modeled  with  different  p-doping  conditions,  as  shown  in  Fig.  4.  We  use  a  uniform  doping  in 
the  absorber,  which  can  be  well  controlled  during  crystal  growth.  When  the  absorber  region  is 
lightly  p-doped,  the  quantum  efficiency  decreases  slightly  because  more  holes  are  introduced 
in  the  VB.  But  the  PL  intensity  increases  with  p-doping  since  the  VB  becomes  increasingly 
occupied  by  holes  (the  quasi  Fermi-level  for  holes  shifts  to  lower  energy),  and  allows  a  larger 
spontaneous  emission  rate.  Since  the  doping  is  light  and  uniform,  the  problem  of  band  distortion 
should  not  be  significant.  The  fact  that  the  light  p-doping  does  not  change  the  cutoff  wavelength 
and  only  slightly  affects  the  QE,  one  can  use  a  lightly  doped  absorber  region  to  reduce  the  dark 
current  and  drastically  improve  the  electrical  performance  of  photodetectors. 

Our  theoretical  results  are  compared  with  the  experimental  data  in  Fig.  5  and  Fig.  6,  show¬ 
ing  a  very  good  agreement.  The  measured  absorption  spectra  are  converted  from  the  quantum 
efficiency  (QE)  spectra  [19,20]  for  a  mid-wave  infrared  (MWIR)  detector  with  a  3.6  pm  cutoff 
wavelength  and  a  long-wave  infrared  detector  (LWIR)  with  a  12  pm  cutoff  wavelength.  The 
absorption  coefficient  is  an  intrinsic  property  of  the  designed  T2SL  while  quantum  efficiency  is 
a  device-dependent  property,  which  can  be  affected  by  doping  concentration,  device  thickness, 
minority  carrier  lifetime,  minority  carrier  diffusion  length,  and  so  on.  The  absorption  coeffi¬ 
cient  can  be  measured  from  the  absorption  difference  between  two  identical  T2SLs  except  with 


for  p-type, 

for  n-type, 
for  intrinsic. 


(22) 
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Fig.  3.  Theoretical  TE  absorption  spectra  of  an  InAs/GaSb  44  A/21  A  T2SL  on  GaSb 
substrate  using  Dirichlet  boundary  condition  (DBC)  with  7  (blue  solid),  9  (green  solid), 
and  1 1  (red  solid)  periods,  and  using  the  periodic  boundary  condition  (PBC)  with  1 1  (red 
dashed),  15  (green  dashed),  and  19  (blue  dashed)  periods.  Inset:  Total  absorption  (blue), 
first  heavy-hole  (HH1)  to  all  conduction  subbands  (CBs)  absorption  (black),  first  light-hole 
(LH1)  to  all  CBs  absorption  (green),  and  second  heavy-hole  (HH2)  to  all  CBs  absorption 
(red). 


Fig.  4.  Modeled  external  quantum  efficiency  and  photoluminescence  (PL)  spectra  of  a  300- 
period  lightly  p-doped  InAs/GaSb  45  A/24  A  T2SL  absorber  with  NA  =  5  x  1016  cm”3 
(cross),  NA  =  2x  1016  cm”3  (dashed),  NA  =  5  x  1015  cm”3  (solid). 


different  thicknesses,  and  the  calibration  is  automatically  included  during  the  spectral  response 
subtraction.  The  minor  spectral  dips  in  the  measured  QE  spectra  in  Fig.  5  and  Fig.  6  are  due  to 
the  atmospheric  absorption.  The  absorption  data  converted  from  measured  QE  follows  our  the¬ 
oretical  prediction  for  both  MWIR  and  LWIR  cases.  The  error  becomes  larger  when  wavelength 
approaches  cutoff  partially  due  to  the  increasing  difficulty  in  separating  weak  signal  from  large 
noise  fluctuation  during  measurement. 

Small  variations  in  material  parameters  and  the  layer  structure  become  critical  when  model¬ 
ing  optical  properties  at  long  wavelength.  A  slight  fluctuation  in  the  input  parameters,  such  as 
the  layer  thicknesses  and  the  valence  band  offset  (VBO),  can  cause  a  large  percentage  change 
with  respect  to  the  small  transition  energy.  In  LWIR  range,  even  5  meV  change  in  the  VBO 
results  in  nearly  1  /i m  shift  in  cutoff  wavelength,  as  shown  in  Fig.  6. 
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Fig.  5.  Comparison  between  theoretical  model  and  experimental  data  [19]  for  an  InAs/GaSb 
18  A/22  A  MWIR  T2SL  on  GaSb  substrate.  Both  the  periodic  boundary  condition  (PBC) 
and  Dirichlet  boundary  condition  (DBC)  (inset)  are  used,  along  with  different  linewidth 
parameter  y. 


Wavelength  (urn) 

Fig.  6.  Comparison  between  our  theoretical  model  and  experimental  data  [20]  for  an 
InAs/GaSb  48  A/22  A  LWIR  T2SL  on  GaSb  substrate.  Both  the  periodic  boundary  condi¬ 
tion  (PBC)  and  Dirichlet  boundary  condition  (DBC)  are  used  and  different  valence  band 
offset  (VBO)  values  take  into  account  the  interfacial  effect  if  actual  interfacial  layers  are 
not  included  in  the  model. 


One  major  difficulty  of  modeling  semiconductor  nanostructures  with  heterojunctions  is  the 
determination  of  the  band  alignment  and  the  VBO.  One  commonly  used  method  is  the  model- 
solid  theory  [30],  where  the  unstrained  valence  band  edges  for  all  materials  are  measured  from 
an  absolute  reference.  However,  a  large  variation  exists  among  reported  experimental  VBO 
values  obtained  through  different  measuring  techniques  or  by  different  researchers.  Besides 
experimental  error,  one  of  the  major  reasons  is  the  imperfectness  of  the  heterojunction  formed 
in  practice.  Interfacial  layers  are  formed  unintentionally  at  the  heterojunction  and  will  affect 
the  measurement  of  VBO  and  other  optical  properties. 

In  this  paper,  we  include  actual  IF  layers  in  our  model  to  study  their  effect  on  optical  prop¬ 
erties.  IF  layers  are  forced  at  the  heterojunction  while  the  thickness  of  one  superlattice  period 
remains  the  same.  If  the  interfacial  type  between  InAs  and  GaSb  is  InSb,  the  unstrained  CB 
edge  for  InSb  is  in  between  that  of  InAs  and  GaSb,  and  the  unstrained  VB  edge  is  slightly 
higher  than  that  of  GaSb,  as  shown  in  Fig.  7  inset.  However,  due  to  the  large  lattice  mismatch 
(~  6%)  between  InSb  layer  and  GaSb  substrate,  the  strain  effect  plays  an  important  role  in  the 
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Fig.  7.  Modeled  TE  absorption  (cross)  and  PL  (solid)  spectra  for  InAs/GaSb  44  A/21  A 
T2SL  on  GaSb  substrate  when  including  InSb  interfacial  layers  of  different  thicknesses:  0 
A  (purple),  1  A  (red),  1.5  A  (green),  2  A  (blue),  2.5  A  (black).  The  inset  shows  the  band 
edge  profile  at  InAs/GaSb  interface  with  forced  InSb  layer  in  between.  The  blue  and  black 
lines  are  unstrained  conduction  band  and  valence  band  edges  respectively.  The  green  and 
red  dashed  lines  are  conduction  and  heavy-hole  band  edges  respectively  including  strain 
effect. 


band  alignment.  The  large  biaxial  compressive  strain  shifts  the  CB  edge  in  InSb  higher  by  377 
meV,  and  also  shifts  the  HH  edge  higher  by  230  meV. 

Assuming  pure  InSb  IF  layers  can  be  intentionally  grown  between  the  InAs  and  GaSb  layers, 
the  VB  states  in  the  T2SL  are  shifted  to  higher  energies  since  the  VB  edge  in  InSb  is  much 
higher  than  the  VB  edges  in  both  InAs  and  GaSb.  The  change  in  the  CB  state  energies  is  small 
compared  to  the  change  in  the  VB  state  energies  because  the  CB  edge  in  InSb  is  still  between 
the  CB  edges  in  InAs  and  GaSb.  As  a  result,  the  InSb  IF  layer  shrinks  the  effective  band  gap  and 
increases  the  transition  wavelength.  Such  effect  can  be  observed  from  the  red-shift  of  the  PL 
peak  and  the  absorption  cutoff  wavelength  in  Fig.  7.  In  our  model,  we  include  IF  layers  at  both 
InAs-on-GaSb  and  GaSb-on-InAs  interfaces.  Practically,  depositing  InSb  on  InAs  interfaces  is 
much  harder  than  that  on  GaSb.  If  the  IF  layer  is  only  deposited  on  one  type  of  interfaces,  the 
effect  of  wavelength  shift  becomes  approximately  half  of  that  on  both  interfaces. 

There  is  normally  a  tradeoff  between  the  cutoff  wavelength  and  the  absorption  strength  when 
designing  T2SL,  but  the  InSb  IF  layers  are  able  to  simultaneously  increase  both  properties,  as 
shown  in  Fig.  7.  The  absorption  coefficient  increases  with  InSb  IF  layer  thickness  because  it 
improves  the  wavefunction  overlap  between  electrons  and  holes.  Electrons  tunnel  more  easily 
through  the  interface  and  the  hole-confining  regions  are  effectively  wider.  As  a  result,  the  InSb 
IF  layer  becomes  a  preferable  region  for  electrons  and  holes  to  recombine.  In  reality,  GaAs 
or  even  ternary  and  quaternary  compound  IF  layers  can  also  be  forced  intentionally  or  formed 
unintentionally  during  the  crystal  growth  of  the  InAs/GaSb  T2SL.  Due  to  the  large  biaxial 
tensile  strain  in  GaAs  on  GaSb  substrate,  optical  properties  are  affected  in  the  opposite  way  as 
by  InSb  IF  layers.  A  wide  tunable  range  of  optical  properties  can  be  achieved  by  controlling 
both  the  composition  and  the  thickness  of  the  IF  layers. 

4.  Conclusion 

We  have  presented  a  comprehensive  model  of  the  electronic  band  structure  and  optical  proper¬ 
ties  of  T2SL  photodetectors  using  the  8-band  k  ■  p  method.  The  theoretical  formulation  for  the 
8-band  Hamiltonian  and  the  explicit  expressions  of  the  8-band  momentum  matrix  elements  are 
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included.  Material  parameters  of  Sb-based  binary  compounds  commonly  used  in  type-II  struc¬ 
tures  are  updated  and  summarized.  Numerical  results  using  two  different  boundary  conditions 
are  compared,  showing  good  agreement.  Our  theoretical  results  of  the  T2SL  optical  properties, 
including  the  absorption  and  photoluminescence  spectra,  agree  very  well  with  the  experimen¬ 
tal  data.  Our  model  also  explains  well  the  doping-dependence  of  the  quantum  efficiency  and 
photoluminescence  in  T2SL  photodetectors. 

The  interfacial  (IF)  effect  in  T2SLs  is  investigated  using  our  model.  Instead  of  adjusting  VBO 
to  account  for  the  IF  effect,  we  include  the  actual  IF  layers  in  our  model.  IF  layers  are  shown  to 
have  a  large  influence  on  the  optical  properties  in  T2SLs.  It  is  predicted  that  by  precise  control 
of  the  IF  layer  component  and  thickness,  we  can  achieve  a  large  tunable  range  of  the  optical 
properties  in  T2SLs  and  significant  improvement  of  the  T2SL  device  performance. 

A.  The  8x8  Luttinger-Kohn  Hamiltonian  in  the  k  p  method 

After  the  basis  transformation  and  axial  approximation,  the  upper  and  lower  parts  of  the  block- 
diagonalized  Hamiltonian  [9]  can  be  written  as, 
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and  m'c ,  71,  72,  73  are  the  corrected  electron  effective  mass  and  corrected  Luttinger  parameters 
based  on  experimental  values  and  Eq.  (7)-(8).  Pcv  is  the  Kane’s  parameter  corresponding  to  the 
momentum  matrix  element  for  s-state  to  p-state  interband  transitions,  which  is  also  represented 
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by  the  energy  parameter  Ep.  ac,  av,  and  b  are  the  Pikus-Bir  deformation  potentials  that  introduce 
strain  effect  into  the  Hamiltonian.  For  quantum  well  structures,  kz  is  replaced  by  the  differential 
operator  —id/dz  in  the  Hamiltonian. 

The  set  of  bases  which  block-diagonalizes  the  Hamiltonian  is  chosen  as, 

M  =  -4[l1)«"^/2  +  il5)^/2].  I“s>  = 


V2l 


V2l 
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The  original  set  of  bases  for  the  8x8  LK  Hamiltonian  is: 
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B.  The  8x8  Momentum  Matrix  Elements 


The  8x8  momentum  matrix  [14,  15]  p  can  be  derived  from  the  8x8  Hamiltonian  under 
the  same  set  of  bases,  following  Eq.  (15).  Each  energy  eigen-state  can  be  expressed  in  terms 
of  basis  functions  weighted  by  envelope  functions  using  k  •  p  perturbation.  The  momentum 
matrix  elements  for  the  transitions  between  two  energy  eigen-states  are  then  evaluated.  Our 
contribution  is  to  derive  the  general  and  explicit  momentum  matrix  elements  in  terms  of  the 
envelope  functions,  which  can  be  readily  used  after  solving  the  coupled  Hamiltonian.  The 
results  are  shown  as  follows,  where  (j)  =  tan ~l(ky/kx). 
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(B-l) 

(T,"’°'1_l|x  ■  p|T,,"’ff2~i)  can  be  obtained  by  replacing  g(j„  with  gi+4,n,  gi,m  with  gi+4.m,  and 
d/dz  with  (— d/dz ),  based  on  Eq.  (B-l). 


#156797  -  $15.00  USD  Received  19  Oct  2011;  revised  29  Nov  2011;  accepted  19  Dec  2011;  published  18  Jan  2012 
(C)  2012  OSA  16  January  2012  /  Vol.  20,  No.  2  /  OPTICS  EXPRESS  2333 


(VF"’0'1  ^-pir,"’02  l)  =  i sin <0  J dz\hk,[  ~  V3fglnglm-V6fglng^m  +  V3fg*3ng6„ 

^^Y&4,n&6,m\  8\,n87,m  ^^8l:n8s,m  ^/^82,n85,m  &3,n&5,m 

,  pz  *  I  .  y/^Y3^  r  *  ^  d  *  *  d  d  *  , 

i  V  4 S4,n85,m\  '  ^  Lo2,«  ^^67,ra  87, m  ~^_82,n  '  83, n  ^^&6,m  <^6,m 


\/3y3/z  d  *  *  d  d  d 

'  2^/2  ^’m  dz^4'n  dz°^’m  //7°2,b"'62,«  ^68, mJ 


W 


"<fec 

J 

7c/z^ 


+ 


3^/i 

2v/2^8’mrfz 


*  *  ^  rf  *  ,  *  ^  i\ 

<?3,n  <?3,;i  8l  ,m~^,84,n  '  84,n~^8l ,m\  j  ■ 


(B-2) 

(*P"’°l_i|x- p|'P(T’02_£/)  can  be  obtained  by  replacing  with  g/+4,„,  g;.m  with  gi~4.m,  and 
d/dz  with  (— d/dz ),  based  on  Eq.  (B-2). 
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{’F"’cl!  ,J|z-p|'F"'!°2=/')  can  be  obtained  by  replacing  gLn  with  g,+4,«,  gi,m  with  gl+4.m,  but 
not  replacing  d/dz  with  (—d /dz),  based  on  Eq.  (B-3). 


The  y-polarized  momentum  matrix  elements  are  different  from  the  x-polarized  ones  only  in 
the  ^-dependent  terms.  However,  after  the  0-integration  in  Eq.  (13),  both  x-  and  y-  polarizations 
give  the  same  TE  transition  rate. 
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replacing  sin0  with  (—cos  0)  . 

(xF"’°1~L|y  •  p|xF'"'°2_L)  can  be  obtained  from  the  expression  of  (vF"’°1~z'|x  •  p|vF'"’°2~L)  by 
replacing  cos  0  with  sin  0  . 
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